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Abstract 

At high magnetic field, the semiclassical approximation which underhes the 
Ginzburg-Landau (GL) theory of the mixed state of type II superconductors 
breaks down. In a quasi-lD superconductor (weakly coupled chains system) 
with an open Fermi surface, a high magnetic field stabilizes a cascade of super- 
conducting phases which ends in a strong reentrance of the superconducting 
phase. The superconducting state evolves from a triangular Abrikosov vortex 
lattice in the semiclassical regime towards a Josephson vortex lattice in the 
reentrant phase. We study the properties of these superconducting phases 
from a microscopic model in the mean-field approximation. The critical tem- 
perature is calculated in the quantum limit approximation (QLA) where only 
Cooper logarithmic singularities are retained while less divergent terms are 
ignored. The effects of Pauli pair breaking (PPB) and impurity scattering are 
taken into account. The Gor'kov equations are solved in the same approxi- 
mation but ignoring the PPB effect. We derive the GL expansion of the free 
energy and obtain the specific heat jump at the transition. We find that each 
phase is first paramagnetic and then diamagnetic for increasing field, except 
the reentrant phase which is always paramagnetic. We also show that a gap 
opens at the Fermi level in the quasi-particle excitation spectrum. The QLA 
clearly shows how the system evolves from a quasi-2D and BCS-like behav- 
ior in the reentrant phase towards a gapless behavior at weaker field. The 
calculation is extended beyond the QLA taking into account all the pairing 
channels and the validity of the QLA is discussed in detail. We show that the 
complete excitation spectrum exhibits gaps at, below, and above the Fermi 
level. We also calculate the current distribution. 

PACS numbers: 74.20-z, 74.70-Kn, 74.90+n, 74.60-w 
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I. INTRODUCTION 



The equilibrium state of type II superconductors was first described by Abrikosov using 
a phenomenological Ginzburg-Landau (GL) theory,EI which was later justified by Gor'kov 
in a microscopic model.! The Ginzburg-Landau- Abrikosov- Gor'kov (GLAG) theory treats 
the magnetic field semiclassically and therefore can be justified in clean materials only at 
high temperature or low magnetic field.! In the last few years, there has been a lot of work 
devoted to the theoretical understanding of the effect of magnetic fields on the mean-field 
theory of the superconducting instability from a completely quantum point of view. 

Most of these works have been concerned with the effects of Landau level quantization in 
superconductors with an isotropic dispersion law.iHll On the one hand, the quantum effects 
of the field have been studied in the vicinity of the semiclassical critical field Hc2{T = 0), 
with emphasis on the precise vortex lattice structure, the quasi-particle excitation spectrum, 
and the de Haas-van Alphen oscillations arising in the mixed state as a consequence of 
Landau level quantization. The first observation of these quantum magnetic oscillations in 
the mixed state occurred nearly twenty years ago in the layer compound 2H-NbSe2,0 and 
interest has been renewed recently with their observation in several other materials. On 
the other hand, it has been proposed that Landau level quantization can lead to reentrant 
behavior at very high magnetic field when the cyclotron energy becomes larger than the 
Fermi energy (cUc ^ -Ei?).!'! This effect is absent from the GLAG theory which predicts a 
complete disappearance of the superconducting phase due to the orbital frustration of the 
order parameter in the magnetic field. This reentrant behavior originates in the suppression 
of the orbital frustration when the electrons reside in only one or the few lowest Landau 
levels. Indeed, when only one Landau level is occupied, the supercurrents can be made to 
coincide with the orbital motion of the electrons in this Landau level if the periodicity of the 
vortex lattice is approximately equal to the orbit radius of the lowest Landau level. Moreover, 
in this very high field limit, it has been argued that the destruction of superconductivity by 
the Pauli pair breaking (PPB) effect can be avoided because the effective ID dispersion law 
allows one to construct a Larkin-Ovchinnikov-Fulde-Ferrell (LOFF)0 state which can exist 
far above the Pauli limited field. The reality of this reentrant superconductivity remains 
however controversialllllll and there has been no experimental result up to now. 

The quantum effects of the magnetic field were also studied in the case of quasi- 
one-dimensional superconductors (weakly coupled chains systems) with an open Fermi 
surface^^^ These effects are especially pronounced when the zero-field critical temper- 
ature Tco is smaller (but not much smaller) than the interchain coupling in the direction 
perpendicular to the field (we will only consider this limit in this paper). (The chains are 
parallel to the x axis. The external field is along the y direction and the interchain hopping 
ty in this direction is assumed to be much larger than t^). In this case, the superconductiv- 
ity is well described semiclassically by the anisotropic GL theory. In particular, there is no 
Josephson coupling between chains even at T = 0. Because of the quasi- ID structure of the 
Fermi surface, the semiclassical orbits in the presence of the field are open. Consequently, 
there is no Landau level quantization but the field induces a 3D/2D crossover :l3E§E3 the 
electronic motion remains extended along the chains and along the direction parallel to the 
field, but becomes confined in the z direction with an extension ~ ctz/uJc (c is the interchain 
spacing in the z direction and Uc is the frequency of the semiclassical orbits) . This dimen- 
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sional crossover is at the origin of a very unusual phase diagram. In particular, it leads to a 
restoration of time-reversal symmetry (as far as the Zeeman splitting is ignored) in very high 
magnetic field {uj^ ^ tz) which results in a reentrant behavior of the superconducting phase 
with a Josephson coupling between chains. This high-field-superconductivity can survive 
even in the presence of PPB because the quasi-lD Fermi surface allows one to construct a 
LOFF state (for any value of the magnetic field) which can exist far above the Pauli limited 
field. Although the origin of the reentrant behavior is very different in the quasi-lD case 
and in the isotropic case, in both cases it appears as a consequence of a reduction of dimen- 
sionality, from 3D to 2D in the quasi- ID case, and from 3D to ID in the isotropic case. The 
suppression of the orbital frustration originates in this reduction of dimensionality. 

Besides the qualitative differences between isotropic and quasi-lD superconductors, there 
is also an important quantitative difference. In the isotropic case, the temperature and mag- 
netic ranges where quantum effects are expected to be important are determined by the Fermi 
energy Ep. For this reason, superconductivity is destroyed for intermediate fields, i.e. for 
fields much larger than in the semiclassical regime but much smaller than in the reentrant 
regime. Moreover, the reentrant behavior can be observed only at very low temperatures 
and very high fields. This restricts considerably the possible candidates to the experimen- 
tal observation of very high-field superconductivity and is one of the reasons which explain 
the absence of experimental results. In the quasi-lD case, it is the coupling tz between 
chains which plays the crucial role. Since can be smaller than 10 K in organic conduc- 
tors, the temperature and magnetic field ranges where very high-field superconductivity is 
expected can be experimentally accessible if the appropriate (i.e., sufficiently anisotropic) 
materials are chosen.@'0 This also means that superconductivity can survive even for in- 
termediate fields between the GL and the very high field regimes. The interest in quasi-lD 
superconductors has been recently raised by experimental results on the organic compound 
(TMTSF)2C104.il Resistive measurements have shown an anomalous behavior of the crit- 
ical field Hc2- Although they do not give a definite answer for the existence of high-field 
superconductivity in (TMTSF)2C10a, these results might be interpreted as the signature of 
a high-field superconducting phase.iil'@ 

The main features of the phase diagram of a quasi- ID superconductor are now well 
understood.0^ Between the GL regime (where the superconducting state is a triangular 
Abrikosov vortex lattice) and the reentrant phase (where the superconducting state is a tri- 
angular Josephson vortex lattice), the magnetic field stabilizes a cascade of superconducting 
phases separated by first order transitions. In these quantum phases, the behavior of the 
system (and in particular the periodicity of the order parameter) is not determined any 
more by the (semiclassical) GL coherence length iz{T) ~ but by the transverse (i.e., 

perpendicular to the chains) magnetic length ctz/ojc ^ 1/H {H being the external magnetic 
field). When entering the quantum regime, Uc ~ T, the transverse magnetic length is much 
larger than the GL coherence length ^^(T). This results in an increase of the transverse 
periodicity of the order parameter and in a strong modification of the vortex lattice:^ in 
the quantum regime, the amplitude of the order parameter and the current distribution 
show a symmetry of a laminar type while the vortices still describe a triangular lattice. The 
existence of this somehow new superconducting state is due to the symmetry of the one- 
particle wave-functions which is incompatible with the symmetry of the Abrikosov vortex 
lattice. The cascade of first order phase transitions originates in commensurability effects 
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between the periodicity of the order parameter (i.e., the transverse magnetic length) and 
the crystalhne lattice spacing. 

In this paper, we study the transition line and the properties of the superconducting 
phases in the mean-field approximation startingfrom a microscopic model. Some of the 
results presented here were published elsewhere." We assume that the superconductivity 
is due to an effective attractive electron-electron interaction of the BCS type. We also 
assume that the quasi-lD conductor is well described above the transition line by the Fermi 
liquid theory, which justifies the use of a mean-field theory. This situation will be realized 
if the system undergoes a single particle dimensionality crossover at a temperature Tr^i > 
TcQ. Below Tj.1, the particle-particle (Cooper) and particle-hole (Peierls) channels decouple 
so that the usual mean-field (or ladder) approximation is justified provided that the bare 
parameters of the Hamiltonian are replaced by renormalized ones in order to take into 
account the effects of ID fluctuation^^ (see also Refs. for a discussion of the validity 
of the mean-field approximation, in particular for the organic conductors of the Bechgaard 
salts family). As pointed out by Yakovenko,lli it is very important that the electrons in 
the {x,y) planes have a 2D behavior below T^i. Even if the magnetic field suppresses the 
electron hopping in the z direction, it has no effect on the electron motion in the (x, y) 
planes and the Cooper and Peierls channels remain decoupled.il This does not exclude the 
existence of thermodynamical fluctuations, in particular in the reentrant phase [uc 3> t^) 
where the system becomes effectively quasi-2D. In this very high field limit, the transition 
from the metallic phase towards a phase with real superconducting long-range order might 
be replaced by a Kosterlitz-Thouless transition.lli This aspect will however not be considered 
any further and we will restrict ourselves to the mean-field analysis. 

In the next section, we calculate the eigenstates and the Green's functions of the normal 
phase in the presence of a uniform magnetic field H(0, H, 0). We use the gauge A{Hz, 0, 0) 
which presents the advantage to yield a very clear physical picture of the dimensional 
crossover induced by the magnetic field. In Sec. [m] , we derive the transition line in the quan- 
tum limit approximation (QLA) where only Cooper logarithmic singularities are retained 
while less divergent terms are ignored. Although this approximation strongly underestimates 
the critical temperature, it provides a clear physical picture of the pairing mechanism re- 
sponsible for the superconducting instability. Moreover, the effects of disorder and PPB can 
easily be incorporated in this approach. In Sec. |V|, we study the superconducting phases 
in the QLA ignoring the PPB effect. We first construct a variational order parameter using 
the results of Sec. |T| and then solve the Gor'kov equations. We derive the GL expansion 
of the free energy and obtain the specific heat jump at the transition. The discontinuity of 
the specific heat jump at the first order transitions is related to the slope of the first order 
transition lines. We find that each phase is first paramagnetic and then diamagnetic for 
increasing field, except the reentrant phase will is always paramagnetic.il We also show that 
a gap opens at the Fermi level in the quasi-particle excitation spectrum. The QLA clearly 
shows how the system evolves from a quasi-2D and BCS-like behavior in the reentrant phase 
towards a gapless behavior at weaker field. In Sec. |V|, we go beyond the QLA. We first 
obtain the transition line and construct a variational order parameter, thus recovering the 
results obtained in Ref.t3 in the gauge A'(0, 0, — ifx). We then derive the GL expansion 
of the free energy. We discuss the importance of the screening of the external field by the 
supercurrents and also compare the results with those obtained in the QLA. The quasi- 
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particle excitation spectrum is obtained from the Bogoliubov-de Gennes equations. Besides 
gaps which open at the Fermi level as obtained in the QLA, gaps open below and above the 
Fermi level. This excitation spectrum is very reminiscent of the one of the field-induced-spin- 
density-wave (FISDW) phases which appear when the effective electron-electron interaction 
is repulsive. Finally, the current distribution is calculated. 



II. GREEN'S FUNCTIONS OF THE NORMAL PHASE 



In this section, we derive the Green's functions in the normal metallic phase in the 
presence of a uniform magnetic field H(0, if, 0). Contrary to what can be found in general 
in the literature concerning quasi-lD conductors in a magnetic field, we work in the gauge 
A{Hz,0,0). The one-particle Hamiltonian is obtained from the Peierls substitution Hq — 
E^k — zV — eA). The dispersion law is given hj {h = ks = 1 throughout the paper and 
the Fermi energy Ep is chosen as the origin of the energies) 

E(k) = v{\kx:\ — kp) + ty cos{kyb) + cos{kzc) , (2.1) 

where v is the Fermi velocity for the motion along the chains {x axis) and ty, t-, are the 
couplings between chains separated by the distance h and c. The condition ty^tz <C Ep 
ensures that the Fermi surface is open. Except in a few cases which will be pointed out 
when necessary, we will not explicitly consider the y direction parallel to the magnetic field 
which does not play any role for a linearized dispersion law (as long as Cooper pairs are 
formed with states of opposite momenta in this direction). In order to take into account the 
y direction, we just have to replace the 2D density of states per spin iV(0) = l/7rt>c by its 
3D value l/nvbc. It should be noted here that no generality is lost at the mean-field level 
when studying a 2D system instead of a 3D system. This is due to the fact that the kinetic 
energy mainly comes from the motion along the chains, which is not affected by the field (see 
below), so that the electron-electron interaction can still be treated in perturbation for a 2D 
system. This should be contrasted with 2D isotropic systems in high magnetic field where 
the perturbative treatment (i.e. the mean-field analysis) of the superconducting instability 
is highly questionable. i'El 

Since the magnetic field does not couple the two sheets of the Fermi surface,ii we can 
write an Hamiltonian for each sheet of the Fermi surface: 



Ti-Q^^ = v{—iadx — kp) + amujc + cos{—icdz) + crh , (2.2) 

where a = + (— ) labels the right (left) sheet of the Fermi surface, m is the (discrete) 
position operator in the z direction and a = + (— ) for | (|) spin, ah = ajjisH is the 
Zeeman energy for a g factor equal to 2. We have introduced the energy ujc = Gv where 
G = —eHc is a magnetic wave vector. The operator —idx commutes with the Hamiltonian 
so that the momentum fc^. along the x axis is a good quantum number. We therefore look 
for a solution (j)'^^{x,m) = e^^'^^ (j)1^{m) . The Fourier transform ^^^(/c^) of (f)1^{m) is solution 
of the Schrodinger equation (setting c = 1 for simplicity) 

[iau^dk, + cos(fc^)]0^^ (/c J = e0^^ (fc^) , (2.3) 
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where e = e — v{akx — kp) — crh and e is the eigenenergy of the eigenstate 0^^(a;, m). The 
solution of (12.31) is 



~e-'-[^''^^-*^""('=^)l (2.4) 
up to a normahzation factor. Going back to real space, we obtain 

<j>l (x, m) ~ ^e'^^^^'-^c'^^'^J^ ^'^^'^^ . (2.5) 
^ J —77 27r 

(x, m) is non zero only if e = aluc where / is an integer. Therefore, the normalized 
eigenstates and the eigenenergies of the Hamiltonian Ti.Q^ are given by (writing explicitly 
the interchain spacing c) 

^l(r) = ^e^'^^^^Ji.m{ai), (2.6) 

(^k^,i,a = v{ak^ - kp) + aluc + ah , (2.7) 

where r = (x, m) and J/ is the /th order Bessel function. is the length of the system in 
the X direction and t = t^/ujc is a reduced interchain coupling. The state ; is localized 
around the /th chain with a spatial extension in the z direction of the order of tc which 
corresponds to the amplitude of the semiclassical orbits.^ The advantage of working with 
the vector potential A{Hz, 0, 0) is now clear: in this gauge, the eigenstates are localized, 
which corresponds to the real physical behavior of the particles as can be seen by examining 
the one-particle Green's function in real space.!! Moreover, since the momentum k^ remains 
a good quantum number, the motion along the chains is in some sense trivial. Note that 
the states ; can be deduced by the appropriate gauge transformation from the localized 
states introduced by YakovenkoH in the gauge A'(0, 0, —Hx). 

In the Hamiltonian (|2.2| ), the magnetic field appears only through the additional term 
aTTiuJc- The effect of the magnetic field is therefore similar to an "electric field" —avH 
whose sign would depend on the sheet of the Fermi surface. This fictitious electric field 
introduces an additional "potential energy" amUc on each chain m which competes with the 
hopping term tzCos{—icdz) and tends to prevent the electronic motion in the z direction. 
The localization in the z direction is almost complete when the difference of "potential 
energy" Uc between two chains is much larger than the transfer integral between chains 
{u!c ^ tz). The fact that an electric field can localize the wave functions (as obtained here in 
a tight-binding model) has been known for a long time.0 This effect has recently attracted 
a lot of attention in connection with the studies of semiconductor superlattices submitted 
to an electric field.H The quantized spectrum which results from this localization is known 
as a Wannier-Stark ladder. The semiclassical picture of this effect yields the well-known 
Bloch oscillations of a band electron in an electric field. Continuing this analogy between 
electric field and magnetic field in a quasi-lD conductor, we can interpret the semiclassical 
trajectories z = zo + c(tz/uJc) cos{Gx) obtained from Hq as Bloch oscillations of the electrons 
in the magnetic field. Most of the effects which are induced by the magnetic field in a 
quasi- ID conductor can be understood as the consequence of these Bloch oscillations. 

In the next section, it will be useful to use Green's functions in the representation of the 
states Introducing the creation and annihilation operators b'^l i^^, Vj^^ i^^ of a particle of 
spin a in the state we define the Matsubara Green's function G'^{kx, I, ui) by 
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= (^^-C/,J"'' (2-8) 
where u = Un = ■nT{2n + 1) is a Matsubara frequency. 

III. TRANSITION LINE IN THE QUANTUM LIMIT APPROXIMATION 

The exact mean-field critical temperature has been calculated numerically and dis- 
cussed in detail elsewhere.lllllli^ In this section, we calculate Tc in the QLA, where only 
Cooper logarithmic singularities are retained while less divergent terms are ignored. Al- 
though this approximation yields a critical temperature several orders of magnitude smaller 
than the exact mean-field result, it provides a clear physical picture of the pairing mecha- 
nism responsible for the superconducting instability. It will also allow us to give a simple 
description of the properties of the ordered phase below Tc (see Sec. Moreover, the 

effect of disorder can be easily incorporated at this level of approximation. 

The total Hamiltonian is now Ti^ + TYint where the effective attractive electron-electron 
Hamiltonian is described by the BCS model with coupling constant A > 0: 



^int = -^ E /^'rC'nr)C'^(r)^^(r)C(r). (3.1) 



We use the notation / (Pv = cJ2mI dx and a = —a, a = —a. The '?/'"(r)'s are fermionic 
operators for particles moving on the sheet a of the Fermi surface. The interaction is effective 
only between particles whose energies are within Q of the Fermi level. 

We first note that the superconducting instability can be qualitatively understood from 
the Wannier-Stark ladder ( |2.7] ). In zero-field, time- reversal symmetry ensures that -E'i(k) = 
i?j(— k) so that the pairing at zero total momentum presents the usual (Cooper) logarithmic 
singularity ~ ln(27f2/7rT) (7 ~ 1.781 is the exponential of the Euler constant) which results 
in an instability of the metallic state at a finite temperature Tco- A finite magnetic field 
breaks down time- reversal symmetry. Nonetheless, we still have e^^ -|. = ^q^-k^^i^^i for 
Qx = —{h + h)G (if we ignore the Zeeman splitting). Thus, whatever the value of the 
field, some pairing channels present the Cooper singularity if the total momentum along 
the chain q^. is a multiple of G. This results in logarithmic divergences at low temperature 
in the linearized eap equation, which destabilize the metallic state at a temperature T^, 
(0 < Tc < Tco).E3ll3 This reasoning holds also in presence of the Zeeman splitting if we 
consider pairing at total momentum = —{li + l2)G±2h/v for two bars /i, I2 of the 
Wannier-Stark ladder. The shift ±2h/v of the total pairing momentum, which displaces the 
Fermi surfaces of spin | and spin J, relative to each other, partially compensates the PPB 
effect and yields to the formation of a LOFF state.liilli0 

Besides the most singular channels which present the Cooper singularity, there exist less 
singular channels with singularities ~ In |2f2/nc(Jc| (n 7^ 0) for T ^ cjc as will be shown 
below. In this quantum limit {uc ^ T), a natural approximation consists in retaining only 
the most singular channels. This QLA has been used previously in the mean-field theory of 
isotropic superconductors in a high magnetic field.i 
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It is worth pointing out that the same kind of reasoning can explain the appearance 
of the FISDW phases in the presence of a repulsive electron-electron interaction.iil Even 
if we add to the Hamiltonian a second neighbor hopping term t'^ cos{2kzc) (assumed to be 
large enough to suppress any SDW instability in zero magnetic field), the spectrum remains 
a Wannier-Stark ladder (Eq. ( |2.7| )) although the expression of the eigenstates differ from 
(|]6D. Since e^^ ,^ = -eQ^+k,^i2,a for Qx = + (h - l2)G, some electron-hole pairing 
channels present logarithmic singularities if Q^ — '^kp is a multiple of G. At low temperature 
and high enough field, this leads to an instability of the metallic phase with respect to a 
SDW phase at wave vector = 2kp + NO {N integer). Thus, the gauge A{Hz,0,0) 
provides a very natural picture of the quantized nesting mechanismEl which is at the origin 
of the FISDW phases in quasi-lD conductorsS The QLA approximation has also been used 
in this context where it is known as the single gap approximation (SGA). 01111 



A. Without disorder and without PPB 

We first consider the simplest case where both PPB and disorder are neglected. In 
order to obtain the critical temperature, we consider the two-particle vertex function 
-poa (n, r']^, r2) for a pair of particles on opposite sides of the Fermi surface and with 
opposite spins and Matsubara frequencies. F"" is evaluated in the ladder approximation 
shown diagrammatically in Fig. |I|. We first write the two-particle vertex function in the 
representation of the eigenstates 



r»"'(ri,r2ir;,r;,) = — 



CL S S 2: 



(^1, h', I'l, l^) 



xC/i(ri)C-fc.,«2(r2)0fc;,r 



(3.2) 



Here we have used the fact that the total momentum of the Cooper pair along the chains 
is conserved. In the ladder approximation, the vertex F^" (/i, /'j^, /2) is solution of the 
equation 

Tq" (/l, ^i, ^2) = (^1, «; ^2, «|^mt|^'l, a'; ^25 



- ^ (/i, a; /2, al^intl^i', a"; /2) «") 



(3.3) 



where 



x°"(g. + (/i+M 



T 



Y.G-;\KXMGi\q.-KX-^) 



(3.4) 



is the two-particle propagator in the representation of the states <^%^^i- Using the expression 
( p.8| ) for the Green's function G'^{kx,l,uj), we have 



Nio) 



\7rT J \2j \2 AmT J 



(3.5) 



8 



where \I/ is the digamma function. A^(0) = l/vrfc is the density of states per spin at the 
Fermi leveL Note that xilx) = J^aX'^ilx) is the pair susceptibihty at zero magnetic field 
evaluated at the total momentum g^.. The first term on the right-hand side of ( p.3| ) is the 
two-particle matrix element of the electron-electron interaction Hint: 

{h, a; ^2, a I Hint I ^'l, ^2^ '^')^A:i,+fc2a;,fci,+fc^, 

- A / rf^r <Pl^,^ {vY<Pl^,,^ (r)*0^;^,, (r)0|^,, (r) 
= - A4..+fe.,fci,+.^,«'-'^a''i-'^ /J- f e^('^+'-'i-'2)V,,_,,(2tcosx) J,i_,. (2tcosx) . (3.6) 



It is clear from ( |3.3| ) that Cooper logarithmic singularities x"(0) ~ \'n.{2'~)Vt / t:T) appear 
when the intermediate two-particle state corresponds to + 2L"G = (which imposes to 
be a multiple of G) where L" is the center of gravity of the Cooper pair with respect to the 
z direction. Intermediate states with q^ + 2L"G = nG {n ^ 0) lead to weaker logarithmic 
singularities y^ijiG^ ~ In |2fi/?T,cc;c| for ojc ^ T. The QLA, which is expected to be valid 
for ijjc ^ T, restricts the Hilbert space to the intermediate states such that q^ -l- 2L"G = 0. 
Since qx is a constant of motion, the center of gravity L = (/i + /2)/2 = —qx/2G also becomes 
a constant of motion of the Cooper pair in the QLA. Thus, (p.3|) reduces to 



rr('L),L(/, n = -a^m"' + ^x(o) E ^J''"r£'(l'),L(/", , (3.7) 

^ a", I" 

where I = li — I2 and I' = l[ — describe the relative motion of the pair in the z direction. 
qx{L) = -2LG and VJ^' = a'a'''v;,r with 

Viv= —Ji{2icosx)Ji'{2icosx). (3.8) 

JO 271 



Note that VJX' is independent of the center of gravity L of the Cooper pairs. To eliminate the 



dependence on a and a', we write ^q^lL),Li^^^') — — Aa'a'' Tq^(^L)^L{l,l'). Tq^(^L)^L is solution 
of the equation 



r,4L),L(/, n = V,,v + Ax(0) E VirTq.^^L),L{l", n ■ (3-9) 

The preceding matrix equation is solved by introducing the orthogonal transformation 
Ui^i> which diagonalizes the matrix Vi^ir. (U~^VU)i^ii = Si^i>Vi^i. The matrix Tg^(^L)^L = 
U~^Tq^(^L)^LU is diagonal: 

r,.(.),.(M') = VY3^^^, (3.10) 

and we obtain 

^ V 1 - Ax(0)V^.,z. ^ ^ 

The metallic state becomes instable when a pole appears in the two-particle vertex function. 
Using x(0) = N (0) \ii{2jfl / ttT) , we obtain the critical temperature 
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where V/q.'o highest eigenvalue of the matrix V. The critical temperatures are shown 

in Fig. 1^ for the two highest eigenvalues of V. The parameters used in the numerical 
calculations (Tco = 1.5 K and tz = 20 K) are the same as those of Fig. 1 and Figs. 4-10 
of Ref.lli. Fig. 1^ clearly indicates that there are two lines of instability competing with 
each other and leading to a cascade of first order transitions in agreement with the exact 
mean-field calculation of Tc.El The existence of two lines of instability results from the fact 
that Vi^i' = if / and I' do not have the same parity. Diagonalizing the matrix Vi^i' is then 
equivalent to separately diagonalizing the matrices ¥21,21' and V2i+i,2F+i- In the following, 
we label these two lines by Iq = 0,1 so that VJq^/q = max; V2i+«o,2«+/o- Since 2L = li + I2 and 
li — 12 have the same parity, L is integer (half-integer) for /q = (/q = 1) and can be written 
as L = — /0/2 + p with p integer. Correspondingly, we have qx{L) = (/q — 2p)G. It is clear 
that the instability line Iq corresponds to the instability line Q = IqG which was previously 
obtained in another approach where the magnetic Bloch wave vector Q plays the role of a 
pseudo-momentum for the Cooper pairs in the magnetic field.0 

As can be seen from Fig. 0, calculated in the QLA is several orders of magnitude 
below the exact critical temperature except in the reentrant phase: it has been pointed out 
previously that in general the QLA strongly underestimates the critical temperature.^ In 
the QLA, we neglect intermediate pair states with a center of gravity L" ^ L = —qx{L)/2G. 
Since the one-particle states i are localized, \L" — L\ is bounded by ~ tz/ujc- This means 
that the QLA neglects the logarithmic divergences \n.{2VL/ uj^), \\i{VL/uj^, \\i{2VL/tz)- There 
are therefore two cases where the QLA becomes quantitatively correct. Either ujc > 
(which corresponds to the reentrant phase) so that the only logarithmic divergence to be 
considered is the Cooper singularity ln(27fi/7rT). Or the cutoff energy Q is sufficiently 
low for the condition ~ f2 to hold. In a conventional (isotropic) superconductor where 
the attractive electron-electron interaction is due to the electron-phonon coupling, Q is the 
Debye frequency and the condition tUc ~ f2 can never be satisfied for reasonable values of the 
magnetic field. In a quasi-lD superconductor, it has been argued that Q ~ T^i, where T^i 
is the single particle dimensionality crossover temperature below which the system becomes 
3D.ii The reason is that the superconducting instability cannot develop at energies e > T^i 
(or equivalently at length scales < v/T^i) where the behavior of the system is ID. In organic 
superconductors like the Bechgaard salts, T^i can be of the order of 10 — 30 K, so that the 
condition ~ ^ could be realized in particular cases although this remains quite unlikely. 
Moreover, in the weak coupling limit Tc ^ Q, the QLA can never be quantitatively correct 
when entering the quantum regime {uc ~ T). 

From ( p.ll| ), one can see that the superconducting condensation in the channel qx{L), L, Iq 



corresponds to the following spatial dependence for the order parameter 

A,.(L),L,ar) ~ E«'^^,^oCL+i(r)C(i)_,.,i_i(r) • (3.13) 



Noting that the matrices V and U have a range of the order of t (i.e., Vi^ii, Ui^i' are important 
for |/|, \l'\ < i), one can see that A^^(i) ^ has the form of a strip extended in the direction 
of the chains and localized in the perpendicular direction on a length of the order of ct. 
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This is not surprising since Aq^(i) ^ results from pairing between the locahzed states 
The expression ( 3.13| ) of the order parameter at T^. will be used in Sec. [V to construct a 



variational order parameter describing the ordered phase below T^. 



B. Effect of disorder 

We evaluate the effect of disorder on the critical temperature calculated in the QLA. In 
presence of impurity scattering, the pair propagator appearing in the integral eouation for 
the vertex function F""' has to be modified by self-energy and vertex corrections!^ as shown 
diagrammatically in Fig. ^. In the Born approximation, the self-energy is given byll 

E^^^H = -^sgnM, (3.14) 
Zt 

where 1/r = 27t N {0)niV^^ . Ui is the impurity density and Vi is the strength of the electron- 
impurity interaction which is assumed to be local in real space. The elastic scattering time 
T is not affected by the magnetic field because the density of states per spin A^(e, H) = N{0) 
is magnetic field independent. Thus, self-energy corrections can be taken into account by 
the usual replacement u lj = u + sgn(a;)/2r in the expression of the one-particle Green's 
function. Because of the vertex corrections shown in Fig. ^, the pair propagator n2^°^(/i, I2) 
in the pairing channel qx{L), L is determined by the matrix equation 

ns^-(/i,/2) = 5„„«.xSHo) + noxSno) E ^rrnr^(^3,/2), (3.15) 

where uq = 1/2'kN{0)t. The variables k refer to the relative motion of the Cooper pair in 
the z direction. In writing ( [3.15| ), we have used the fact that in the QLA, the only states 
which are allowed satisfy qx{L) = —2LG where L is the center of gravity of the pair with 
respect to the z direction. is defined by 

xm = ^ E GfiK, L + |, cD)Gf (g.(L) - fc., L - |, -Co) . (3.16) 
Introducing the matrix 

n^(/i,/2) = E («ia2)'^ns^°^(/i,/2), (3.17) 

the matrix equation ( p.l5| ) becomes 

n^(/i, h) = xdO) + UoXdO) E ViuiMh, h) , (3.18) 

h 

where ^(^(O) = HaXZi^)- order to obtain the preceding equation, we used the property 
that Vi^ii is non zero only if / and /' have the same parity. The matrix Iljj is diagonalized by 
the transformation U: 



h 



Km-. (3-19) 
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In presence of impurity scattering, the two-particle vertex function is determined in the 
QLA by the equation 

r^:(L),L(^i,^2) = -wr^iT +^ty.t.t. vrx'^T\iM^Z(i),Li.h^h) . (3.20) 

03,^3 04, u 

The dependence on the indices ai is suppressed by writing r^^^^j^/i, = — Aai^a2 r(/i, /2)- 
Using the property that Ilg" (/, V) is non zero only if I and V have the same parity, we obtain 

r,,(L),L(/i,/2) = Vi,,i, + Vi,,i,Uc.{h,U)r,^^L)AU,h) . (3.21) 

The preceding matrix equation is diagonalized by the transformation U: 

From ( p.l9| ) and (|3.22|) , we see that the appearance of a pole in the two-particle vertex 
function corresponds to 

1 - XV,,T Y: -. ... = . (3.23) 

This equation determines the critical temperature in the pairing channel qx{L),L,l. As in 
the pure system, the highest Tc is obtained for / = Iq defined by Vi^^Iq = maxiVi^i. The result 
( ^.23| ) can be expressed as 

Ar(O)lnf^) =Tf El ^7°^ .n^ -^f Ex.(O), (3.24) 

where T^^^ is the critical temperature in presence of disorder and the critical temperature 
of the pure system given by (|3.12|) . Using 

Xc.O=^^, 3.25 

\Ul)\ 

we obtain 

For Tc — Tf^^ <^ Tc, the preceding equation simplifies in 

Tr - Tf'^ vr 



Tc 8rTc 



(l-n„0- (3.27) 



In the preceding equation, Tc is the critical temperature calculated in the QLA without 
impurity scattering. Since Tc is several orders of magnitude below the exact mean-field 
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value, T^^^ will also be much smaller than the exact value. However, a reasonable estimation 
of T^^^ can be obtained using for the exact value instead of the QLA value. From Eq. 
( p.27| ), one can see that Vi^^Iq, which comes from vertex corrections in the pair propagator, 
tends to reduce the difference Tc — T^'^^. According to ( p.27|) , impurity scattering does not 
affect the critical temperature in the reentrant phase since VJo,'o ~^ ^ when Tc —>■ Tco- This 
is a direct consequence of Anderson's theorem which states that the critical temperature 
is independent of a (weak) disorder for a system with time-reversal symmetry.0 Obviously, 
( p.27| ) restricts the observation of high-field superconductivity to clean superconductors with 
a critical temperature not too small. As pointed out in Ref.Ei, this latter condition advan- 
tages materials with a large anisotropy. The consequences of ( p. 27] ) were discussed in detail 
in Ref.lii in the case of the Bechgaard salts. 



C. Effect of PPB 



We evaluate the effect of PPB on the critical temperature calculated in Sec. [ill A| (but 



ignoring the effect of disorder). The equation (|3.3| ) for the two-particle vertex function 
r^^"'(/i,/2;/'i/2) involves the quantity x {Qx + (/" + l2)G) which is given by (^ with the 
replacement avq^ — ^ ctvQx + 2h. Therefore, logarithmic singularities arise through Xq" ilx + 
(/'Z + ZjOG*) each time we have Qx = —2L"G + qo where L" is the center of gravity of the pair in 
an intermediate state and go = ±2/i/f. In the QLA, we retain only the intermediate states 
corresponding to these logarithmic singularities. If h/ijc is not (and not too close to) an 
integer, the center of gravity L of the pair becomes a constant of motion and is related to the 
total momentum by qx{L) = —2LG + go- The equation which determines the two-particle 
vertex function then reduces to 



n = -^v,f + A E v^:^f x""(go)r^:(i), J/", n . (3.28) 

a", I" 

Following the analysis developed in Sec. [Ill A|, we find that a pole appears in the two-particle 



vertex function when 

l-XVijx{qo)=0. (3.29) 

Using 

for h ^ T, we obtain the critical temperature 

where Vi^^ig = max^V;^;. In the reentrant phase ^ t^), T^o so that ttT^q/ i'jh, 

a result which was obtained in Ref.El. The superconducting phase is always stable at T = 0, 



13 



whatever the value of the magnetic field. This divergence of the orbital critical field if ^ ? (T ) 
at low temperature is a characteristic of the LOFF state in a quasi-lD superconductor.ll2lO 
If h/uc is equal (or close) to an integer, the preceding analysis does not hold any more. 
In this case, the center of gravity of the Cooper pair can take the two values —qx/2G±h/uJc. 
This enlarges the number of accessible intermediate states and should lead to an increase of 
the critical temperature with respect to the result ( |3.31|) . It has been proposed by Lebed'Ell 
that this situation could be reached by tilting the field in the (x, y) plane. This lets the 
Zeeman energy h unchanged but modifies the orbital frequency which becomes Uc cos 6 where 
6 is the angle between the field and the y axis. For certain values of 6, h/{ujcCos9) is an 
integer. The effect of the component = if sin 6* is expected to be small (if 6 is not too 
large), because the critical field H^2{^) parallel to the x axis is very large compared to the 
critical fields in the other directions. 



D. Effect of PPB and disorder 



We evaluate the effects of both PPB and disorder on the critical temperature calculated 
in Sec. [Ill A[ Following the analysis developed in the two preceding sections, we find that 
the critical temperature T^^'^^^ is determined by the equation 



A^(0)ln 



1. 
T? 



(3.32) 



where is the critical temperature calculated in the preceding section. For T^f — Tj'*^'^ -C 
T^, we can expand ( p.32|) to first order in uq (or 1/r), which leads to 



iV(o)in(;^) ^TJ2[xUqo)-xUqo)]+uoVi„i,TJ2xliqo) 



where 



Xuj[qo = — sgnl^^j 
vc 



+ 



1 



lioj Hoj — 4/i 



(3.33) 



(3.34) 



7riV(0) 



r 



1 



4cj2 (2iuj - 4/i)2 



Performing the sum over the Matsubara frequencies, we obtain 



) ~ 16tT 



1 



2 ,/l h 



Srh 



7tT 



(3.35) 



+ ^Im^(l + ^A], (3.36) 



where is the first derivative of the digamma function. For h ^ vrT, we can use "^{z) 
Inz — l/2z for \z\ ^1 to get the following expression to lowest order in 1/r: 



T^P,dis rpP 

c c 

C 



'o,'o 1 c 



p2 



(3.37) 
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As in Sec. p!IIB| , a reasonable estimation of the effect of disorder can be obtained from the 
preceding equation using for the exact mean-field value. The effect of vertex corrections 
is weaker than in the absence of PPB (Eq. ( p. 27] )). Even in the reentrant phase where 
Vig^ig 1, the effect of impurities remains important since time- reversal symmetry is broken 
by the Zeeman splitting and Anderson's theorem does not hold. The sensitivity of the LOFF 
state to elastic impurity scattering has been known for a long time.Ei The consequences of 
( p.37| ) were discussed in detail in Ref.lii in the case of the Bechgaard salts. 



IV. ORDERED PHASE IN THE QUANTUM LIMIT APPROXIMATION 

In this section, we study the ordered phase (T < T^) in the QLA where only the pairing 
channels leading to Cooper logarithmic singularities are retained. We only consider the 
orbital effects of the magnetic field, i.e. we put the g factor equal to zero. Using the results 
obtained in Sec. [111 A| , we first construct a variational order parameter with two unknown 
parameters, the amplitude and the periodicity in the z direction. We solve the Gor'kov 
equations in the QLA and obtain the normal and anomalous Green's functions. We then 
deduce the GL expansion of the free energy, the specific heat jump at the transition and the 
quasi-particle excitation spectrum. 

Following the original approach proposed by Abrikosov,BI we construct the order param- 
eter for T < Tc as a linear combination of the functions Aq^(^L)^L^ig{r) (Eq. (|3.13|) ) describing 
the superconducting condensation in the channel qx{L), LJq. Since Tc does not depend on 
L, the most general linear combination can be written as 

A{r) = J2l{L)A,^iL),L,ioir). (4.1) 



As shown in Sec. [Ill A| , 2L must have the parity of /q- Since Ag^(^Q^L^ig{r) is localized in 
the z direction with an extension of ct, a natural choice for the coefficients jiL) is to take 
7(L) 7^ if L = — /o/2 + pN' {p integer) where the unknown integer A^' is expected to be 
of order t. In order to correctly describe the triangular Josephson vortex lattice in the last 
phase (which has periodicity 2c in the z direction) ,t3 we choose 'j{L) = 7p = 1 (i) for p even 
(odd). This leads to (noting N = 2N') 

A,,^(r) = A^f/2,+Zo,Zo7pe^^'°-^^)^V «+,_(t)Jp^_,_,^_(-t) , (4.2) 
i,p 



where the amplitude A is chosen real. Eq. ( |4.2| ) defines a variational order parameter where 
the two unknown parameters A and have to be determined by minimizing the free energy. 
It can be seen that |A(r)| has periodicity = 2ti/NG and = Nc so that the unit cell 
contains two flux quanta 0o for a particle of charge 2e: ifa^a^ = 20o (when a triangular 
lattice is described with a square unit cell, the unit cell contains two flux quanta). In 
Ref.El, the order parameter was constructed by imposing that it describe both the triangular 
Abrikosov vortex lattice in weak field (cUc <^ T) and the triangular Josephson vortex lattice 
in very strong field {uo^ ^ t^). Both approaches lead to the same order parameter when only 
the Cooper singularities are retained. 
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The order parameter is shown in Figs. for the phases N = 2,4,6,8 and 10. The 
exact (at the mean-field level) order parameter corresponding to Figs. |^, |^, and |^, is shown 
in figures 10, 9 and 8 of Ref.lli. The order parameter calculated in the QLA is a good 
approximation of the exact order parameter for the last phases N = 2, 4 and 6. For larger 
values of A^, there appears significant differences between the approximate and exact results. 
For example, in the phase N = 8 at H = 1.7 T, the order parameter obtained in the QLA 
shows zeros which are not present in the exact calculation. 



A. Gor'kov equations 

The starting point of our analysis is the mean-field (or generalized Hartree-Fock) Hamil- 

+ rt-mt 



tonian Hq + with 



n^! = l[d'rY: [A.(r)*C(r)C(r) + H.c] , (4.3) 

where 

A:(r) = AE(#(r)Cnr)). (4.4) 

a 

A|(r) = — Aj(r) is the variational order parameter defined by (|4.2| ). In order to derive the 
thermodynamics and the excitation spectrum in the superconducting phases, it is necessary 
to determine the normal and anomalous Green's functions 

G^(r,r',r) = -(T.C(r,r)C>0)), 

r) = -{T^rJ{r,T)i:p{r,Q)) , (4.5) 

whose Fourier transforms with respect to the imaginary time r are solutions of the Gor'kov 
equations 

(zu - K,.)G^(r, r', u) - A.(r)F#t(r, r', u) = 5(r - r') , 
(-zc.-K,yF#t(r,r',c.) + A:(r)G^(r,r',c.) =0, (4.6) 

with the self-consistency equation 

A:(r) = AT5:F;t(r,r,^). (4.7) 

a,uj 

In the QLA, the supercurrents vanish as will be shown in Sec. 0. The magnetization 
M (which is parallel to the magnetic field by symmetry arguments) and the fiux density 
B = H+AnM are uniform. In the Gor'kov equations, we should therefore in principle replace 
the external magnetic field H by the uniform fiux density B. However, when minimizing the 
free energy with respect to B, it turns out that the approximation B = H is very accurate 
(see Sec. [IV B| ). Moreover, the magnetization vanishing in the middle of each phase (see 
below), this approximation becomes exact at this point. In this section, we consider that 
the fiux density is equal to the external magnetic field. 
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To exploit the translational symmetry of the order parameter, we introduce the magnetic 
Bloch states 0q ; obtained by diagonahzing simultaneously the three operators TYq, e"''^'' and 
which define the magnetic translation group: 



N 



AT ^ 



^q^+pNG,l-pN ■ 



(4.8) 



These states are eigenstates of the Hamiltonian Hq^ with the eigenenergies = e^^ i^^. 
Nz is the number of chains, q is restricted to the magnetic Brillouin zone (MBZ) which is 
chosen to be equal to 



71 71 

'kp — IqG , —kp — IqG H 

dry Qi-y 



for the left sheet of the Fermi surface, and to 



TT TT 



71 , 71 

kp ,kp-\ 



TT TT 



(4.9) 



(4.10) 



for the right sheet of the Fermi surface. This choice is made for later convenience. There 
are N branches crossing the Fermi level on each sheet of the Fermi surface. The excitation 
spectrum of the normal phase in the representation of the magnetic Bloch states is shown in 
Fig. ^ We introduce Green's functions in the representation of these magnetic Bloch states 
writing 



G:(r,r',c.)=y0?(r)0^(r')*G:(l,2, 



1,2 



F;t(r, r', uj)=Y: 0?(r)*0^(rO*F;t(l, 2, 



(4.11) 



1,2 



where we use the notation i = (qj, /j). In this representation, the Gor'kov equations become 
{zoo - 6t)G:(l, 2, ^) - y A^(l, 3)F^^(3, 2, uj) = 6,,, , 

3 

i-zuj - e%)F^\l, 2,uj) + Y: A:(1, 3)*G:(3, 2, ^) = , (4.12) 

3 

where the pairing amplitude A"(l, 2) is defined by 

A^(l,2) = J rf2r</,-(r)>f(r)*A,(r) . (4.13) 

Multiplying ( |4.7| ) by Ao-(r) and summing over r, we can rewrite the self-consistency equation 
as 



/ dh\A^{r)\' = AT^^ A:(l,2)F;t(l,2, 

w,o 1,2 



(4.14) 



A careful examination shows that the pairing between 0^ ^ and 0/J)G-q -lo-i (^i^h G = 
(G, 0)) is the only pairing compatible with the spatial dependence of the order parameter 
([4.21). Thus, in the QLA, we have 
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^"(qi) ^i; q2, h) — '^qi+q2,/oG^«i+i2,-«o^"(qi5 ^i) • 



(4.15) 



The preceding result can also be obtained noting that in the QLA the electron-electron 
interaction conserves the center of gravity of the particles with respect to the z direction 
(Sec. pT]) . Eq. ( |4.15| ) follows from this property of the interaction. The definition of the 
MBZ ( ^ , ^10|) ensures that if q belongs to the MBZ, then ZoG-q also belongs to the MBZ. 
The pairing amplitude A"(q, /) can be calculated using ( pl6| , ^l2| , ^y8| ) : 



A:(q,0 = A^(g„/) 



(4.16) 



and A| = — Aj = A. The equality ([4.15|) and the Gor'kov equations ( [4.12|) imply that the 
Green's functions in the representation of the magnetic Bloch states are diagonal: 



G"{qi, h] q2, k, t^) = '^qi,q2^h,«2<^"(qi> h, uj) , 

^"^qi. ^i; q2, ^) = 5qi+q2,ioG5«i+z2,-io-C^(qi) ^) 

where G"(q, /,cij) and F°"l'(q, /, cu) are given by: 



(4.17) 



-q,/,CT 



u;2 + e°,/ + |A-(q,/)|^ 



A^(q,/)' 



a;2 + e-,/ + |A^(q,0P 



(4.18) 



In the next sections, we will use the preceding expressions of the Green's functions to obtain 
the GL expansion of the free energy and the excitation spectrum. 



B. Ginzburg-Landau expansion and thermodynamics 



We first note that the sign of the magnetization in the ordered phase can be obtained from 
general thermodynamics arguments.^ The magnetization is obtained from M = —dFe/dB, 
where Fe{T,B) is the electronic contribution to the difference of the free energies of the 
normal and superconducting phases. Since Tc is determined by Fe{Tc, B) = 0, the magneti- 
zation M close to the transition line has the sign of dTc/dH. Each phase will therefore first 
be paramagnetic and then diamagnetic for increasing field (which implies that there is a 
value of the field in each phase for which the magnetization vanishes), except the reentrant 
phase which is always paramagnetic.il 

Since the flux density B is uniform in the QLA, the free energy densitv is equal to 
F{T, B) = F(,(T, B) + B'^ /Stx where the electronic contribution can be written^ 



Fe(T,S) = /^rfA'-^A'' 
^ ' ^ Jo dA' 



A.(r) 



A 



2^ A'dA'{X-'-g{A')) J 



r d'^r 


A.(r) 


1 'Y 


A 



(4.19) 
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where the function g{A) = 1/A is defined by ( [4.14|) . In ( [4.19|) , all quantities are calculated 
in presence of a constant fiux density B. In order to get the free energy per volume unit 
(and not per surface unit), the 2D density of states iV(0) = l/nvc should be replaced by its 
3D expression l/nvhc. Expanding in (|4.14| ) in power of A leads to 



(4.20) 



with 



a 



P 



A.(r) 



A 



E 



|A?(q,0/Ap 



T 



T |Ag(q,/)/Ar 



(4.21) 



In the two preceding equations (and in the following), the sum over q is restricted to the 
MBZ. Using ( ^ , gl6|) , and the property 



A^(g.,/ + piV)=e2*^'^^'^^A-(g„/), 



(4.22) 



we obtain 



2V, 



a 



N 
2iV(0) 



N 



Tr 



P = /3- 



BCS 



T/4 N 

7pi7p2 7p37pi+p2-p3 2/+«o+piAf,'o 



1=1 pi,P2,P3 

^U2l+lo+p2N,loU2l+lo+p-sN,loU2l+lo + (pi+P2-P3)N,lo 5 



(4.23) 



where /? is calculated at Tc and Pbcs = 7C(3)A^(0)/87r2T^2. The equilibrium state is deter- 
mined by the minimum of the Gibbs free energy 



G{T, H) = F{T, B) - 
= F,{T,B) + 



BH 

{B - Hf 
8^^ 



Stt 



(4.24) 



Minimizing G{T,H) (or equivalently Ff.{T,B)) with respect to A, we obtain 



a 



(4,25) 



Since a and j3 are known, we can minimize numerically G(T, H) with respect to B and 
then look for the integer N which minimizes the Gibbs free energy. It turns out that 
the minimum is reached for B = H with a very high accuracy. As pointed out above, 
this approximation becomes exact in the middle of each phase where the magnetization 
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vanishes. In the following, we therefore consider that the flux density is equal to the external 
magnetic field. Thus, we can consider the free energy Fg{T,H) instead of the Gibbs free 
energy G{T, H). In the reentrant phase, we find that the minimum of Fe{T, H) is obtained 
for N = 2. When the field is decreased from its value in the reentrant phase, the system 
undergoes a first order transition and the minimum of Fe is then obtained for = 4. This 
is in agreement with Refs.llilli where it is argued that the first order phase transitions are 
due to commensurability effects between the crystalline lattice spacing and the periodicity 
of the order parameter. Unlike what is expected, the best value of N switches to 6 before 
reaching the next first order transition. This indicates that the QLA gives a bad estimate 
of the free energy in the phases > 4. A correct calculation of the free energy in these 
phases requires the inclusion of all the pairing channels and maybe also of the screening of 
the magnetic field. This point will be further discussed in Sec. |V lj| . 

The calculation of the free energy can easily be extended to the case of a more general 
order parameter defined by ( |4.2| ) but with 7p = 1 (7) for p even (odd) where 7 is an unknown 
parameter. For a fixed periodicity N, we have found that the free energy is stationary with 
respect to 7 for 7 = 0, 1 or z which corresponds to square lattices of periodicities = Nc and 
ttz = Nc/2, and to a triangular lattice of periodicity = Nc. Numerically, we have verified 
that in the reentrant phase, the triangular lattice with periodicity N = 2 has a lower free 
energy than the square lattice with periodicity = 1 or iV = 2. This result (stationarity 
of the free energy for the square and triangular lattices, the latter corresponding to the 
minimum of the free energy) extends to the quantum regime Uc ^ T a. property which is 
well-known in the GL regime .i@ 

The specific heat jump at the transition is obtained from AC = —Td'^Fe/dT'^. From 
( ^2^ , ^251) , we obtain 



AC/C_ 



N 



BCS 

N 

^ J2 7pi7p2 7p37pi+p2-p3^2«+/0+Pl^,/0 

1=1P1,P2,P3 

-1 

X U 21+I0 +P2 N,lo U 2I+I0 +P3N,lo U 2I+I0 + (pi +P2 -P3)N,lo 



(4.26) 



where is the specific heat of the normal state and (AC/Civ)Bcs = 12/7^(3) is the zero- 
field value. The ratio r is always smaller than 1 (but close to 1 in the reentrant phase) 
and discontinuous at each first order phase transition (Fig. 0). As in the case of the 
FISDW phases,@ these discontinuities can be related to the slopes AT /AH of the first 
order transition lines. Consider two consecutive superconducting phases denoted by 1 and 
2. For each phase i = 1, 2, the free energy can be written as (T < TJ*)) 

Fi{T, H) = -a,(Tf>, H){T - T^f , (4.27) 

and the specific heat jump at the transition is AQ = 2ai{T^\H)T^\ denotes the 
transition temperature between the normal phase and the superconducting phase i. The 
first order transition line is obtained by writing the equality of the free energy between the 
two phases: 
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- ai{T^^\H){T - TPf = -a2{TP, H){T - T^f 



(4.28) 



Introducing the quantities AT = T — T* and AH = H — H* where the point {H*,T*) is 
defined by the intersection of the two transition fines, i.e. T* = T^^\H*) = T^'^^H*), tfie 
preceding equafity is rewritten as 



'a.iT->, H)[^^-^)^ ^a.iTy>, H)^-^-^j. (4.29) 
Taking tfie fimit T —>■ T* and H H*, we obtain 



AT y/AUT^-y/M^^ 

AH ySUT-v^ ' 

wfiere ACi = ACi{H*) and AC2 = ACiiH*). Since dT^^^/dH < and dT^/dH > 
(Fig. it is clear tfiat tfie slope AT / AH of tfie first order transition line fias tfie sign of 
AC2 — ACi. tfie slope is positive for tfie transition between tfie pfiases = 4 and N = 2 
and negative for tfie otfier transitions. 

Tfie calculation of tfie specific fieat jump at tfie transition clearly sfiows an important 
difference between tfie studv of tfie superconducting pfiases in tfie QLA and tfie study of tfie 
FISDW pfiases in tfie SGA.iii3i In tfie SGA, tfie electron-fiole pairing amplitudes are given 
by a single number (corresponding to tfie gap wfiicfi opens at tfie Fermi level) so tfiat tfie 
tfiermodynamics of tfie FISDW reduces to tfiat of tfie (zero-field) BCS model. For example, 
one obtains r = 1 for any value of tfie magnetic field. In order to get corrections to tfie 
BCS description, one fias to go beyond tfie SGA. In tfie QLA used fiere for tfie study of tfie 
superconducting pfiases, tfie (electron-electron) pairing amplitudes A'^^q^J) do not reduce 
to a single number so tfiat tfie tfiermodynamics is different from tfiat of tfie (zero-field) BCS 
model. 



C. Excitation spectrum 



From tfie Green's functions ( [4.18|) , we deduce tfie quasi-particle spectrum in tfie QLA: 



E. 



± 



± 



(4.31) 



For eacfi brancfi / wfiicfi crosses tfie Fermi level, a gap 2|A°(g2,Z)| opens at tfie Fermi level 
for = akp — IG wfiere ^ g. = 0. Tfie pairing amplitudes A°(g2,Z) are sfiown in Fig. 
[n] for tfie last five pfiases /v' = 2, = 4, iV = 6, = 8 and = 10. Note tfiat 
|A°(— g^,/ ± A^/2)| = |A"(g^,/)|. In tfie very fiigfi field limit t <^ 1, tfie dispersion of 
A"((3'2,/) witfi respect to qz is very weak (of tfie order of t^). Tfie corresponding excitation 
spectrum for tfie rigfit sfieet of tfie Fermi surface is sfiown in Fig. |T2| for tfie two brancfies 
/ = and 1 = 1. Comparing Figs. ^ and |I2|, we clearly see fiow tfie excitation spectrum is 
modified by tfie superconducting pairing. Tfie dispersion of tfie pairing amplitude A"(g^, /) 
witfi respect to q^ increases witfi decreasing field. For large A^ and some particular values of 
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the field (for example in the phase A^ = 8ati^^ = 1.7T), A^(g^, /) vanishes at = ±7r/2a2 
which results in a gapless excitation spectrum. It is rather surprising that the phase N = 8 
is gapless at if = 1.7 T while the phase = 10 is never gapless. This point, which seems 
to indicate again that the validity of the QLA is restricted to the last few phases, will be 
further discussed in Sec. |V C| . Nevertheless, the QLA clearly shows the general tendency: 
the system evolves from a quasi-2D (BCS-like) behavior in very high magnetic field {t ^ 1) 
towards a gapless (or at least with a weak minimum excitation energy) behavior at weaker 
field. 

In the QLA, the zeros of the excitation spectrum are related to the zeros of the order 
parameter in real space. This results from the relation 

(—I)*" 



where = or tt depending on the value of the magnetic field. Eq. ( [4. 321) will be derived 
in Sec. 0. Using (|4.32| ), the order parameter (|4.2|) can be rewritten as 



A,„,^(r) = Ae-*VH„/o(-l)"e^'°''' E lpe-'''''''U,^-pN+io,io ■ (4-33) 



Comparing the preceding equation with (|4.16| ), we see that 



This shows that a gapless spectrum at = ±Ti/2az in the /th branch implies that the 
order parameter in real space vanishes at the points (x = =Fa3;/4,m = I). Such a relation 
between the zeros of the excitation spectrum and the nodes of the order parameter has 
also been obtained by Dukan and Tesanovic in their study of the vortex lattice of isotropic 
superconductors in the QLA.I 

Unlike for thermodynamics quantities, the third direction has a non trivial effect on the 
excitation spectrum. For a 3D system, instead of ( 4.31 ), we have 



Elky,i,. = ±^jK,i,^ + cos(fc,6))2 + I A°(g„ , (4.35) 

where ky is the momentum along the y direction {—i^/h < ky < Tc/b). A'^{qz, I) is the pairing 
amplitude between the magnetic Bloch states 4>q^ky,i,a and 0i(,G-q,-fcy,-io-«,CT- It remains 
unchanged since the Cooper pairs are formed with states of opposite momenta in the y 
direction. The effect of the third direction can be viewed as a shift of the MBZ which is now 
centered around ±kF{ky) defined by 

vikpiky) — kp) + ty cos{kyb) = . (4.36) 

Each time A"(g^,/) vanishes (for a given Bloch momentum and a given branch /), there 
is a line of zeros in the excitation spectrum determined by 

el^i^^ + tycos{kyb)=0, (4.37) 
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where qx belongs to the MBZ and — vr/fe < ky < njh. When calculating thermodynam- 
ics quantities, one has to sum over the quantum numbers associated with the excitation 
spectrum. For a linearized dispersion law (Eq. ( |2.1| )), the ky dependence of the spectrum 
disappears when summing over the momentum and only yields a factor 1/6. This property 
results from the fact that q is coupled only with ZqG — q. Thus, the direction along the field 
modifies the quasi-particle excitation spectrum but not the thermodynamics quantities. 



V. BEYOND THE QUANTUM LIMIT APPROXIMATION 



In the two preceding sections (Sees. |T| and we have studied the quantum super- 
conducting phases in the QLA. This approximation allowed us to obtain a clear and simple 
physical picture. Most of the quantities of interest (critical temperature. Green's functions, 
thermodynamics quantities...) could be expressed in a rather simple form. Nevertheless, 
the QLA suffers from a lot of weak points. It is restricted to the quantum regime Uc ^ T 
and strongly underestimates the critical temperature. It fails to calculate the free energy 
in the phases > 4. Moreover, as will be shown in this section, it does not correctly 
describe the entire excitation spectrum and the supercurrents vanish in this approximation. 
It is therefore necessary to extend the analysis of the two preceding sections to take into 
account all the pairing channels. In the QLA, we only considered the pairing amplitudes 
A"(q, /; IqG — q, —Iq — I). In this section, we will call these pairing amplitudes primary gaps. 
The other pairing amplitudes will be called secondary gaps. The reason for this designation 
is that the primary gaps open at the Fermi level while the secondary gaps open above and 
below the Fermi level (Sec. | V (J| ) . Thus, in this section, we will study how the secondary 
gaps modify the results obtained in the QLA. In the following, we first obtain the exact (at 
the mean-field level) transition line and construct a variational order parameter, thus recov- 
ering the results obtained in Ref.0. Then, from the Gor'kov equations, we derive the GL 
expansion of the free energy, and the excitation spectrum is obtained from the Bogoliubov- 
de Gennes (BdG) equations. Finally, we calculate the current distribution in the ordered 
phase. As in Sec. the PPB effect is ignored. 

A. Transition line and variational order parameter 

In order to obtain the exact (at the mean-field level) transition line, we consider the 
linearized gap equation in a mixed representation {qx,m) by taking the Fourier transform 
with respect to the x direction: 

X'^A{qx,m) = c'^K{qx,m,m')A{qx,m') , (5.1) 



K{qx,m,m') =TJ2 I dx e-'''-^''-'''^G^r,r' ,iu)Gfir,r' , -co) . (5.2) 
Writing the real space Green's function G'^{r, r', u) as 
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G:{t, t', uj) = Y: <PtJr)r,Jr'rG':iK, I, u) , (5.3) 

where (/)^^ ^ and G'^{kx, l.oj) are given in ( |2.6|) and (|2.8| ), we obtain 

cK{qx, m, m') = (-1)—™' ^ Vi^2n.,i~2rn'X{q. + IG) • (5.4) 



The matrix VJ,/' is defined by ( |3.8|) and x(gx) = Z]aX"(Q'x) is the pair susceptibihty at zero 
magnetic field evaluated at total momentum (Eq. ( p.5|) ). It is clear from ( |5.4| ) that the 
highest Tc is obtained when q^ is a multiple of G due to the appearance in this case of Cooper 
logarithmic singularities x(0)- Since the change qx ^ qx ^ 2G is equivalent to the change 
m ^ m ± 1, it is sufficient to consider the two cases = IqG where Iq = 0, 1. Comparing 
(13,0) with Eq. (18) of Ref.S, we note that A(g^ = IqG, m) = Ag^'"^. The coefficients 
were introduced in Ref.lli through the Fourier expansion of the order parameter. Thus, we 
have rederived the linearized gap equation which was obtained in the gauge A'(0, 0, —Hx). 

In the QLA, where only logarithmic singularities are retained, the linearized gap equation 
reduces to 

A-iaQLA(^^ = l^Q^^) = ^(0) Y.{-ir-"''V2m+io,2m'+kA'^''\qx = kG,m') . (5.5) 

m' 

This matrix equation is solved by the orthogonal transformation U which was introduced in 
Sec. m and which diagonalizes the matrix Vi^ii. Obviously, we recover the expression (p.l2|) 
of the critical temperature and 

AQLA(g^ = l^G.m) = {-irU2m+i,M ■ (5.6) 

For qx = {Iq — 2p)G, the order parameter is expressed as 

AQLA(3;, m) = i-ir-^U2m-2p+io,ioe'^''-'''^'''' . (5.7) 

Comparing this expression with ( p.l3| ), we note that we must have 

U2m+io,io = const (-!)"" U2l+lo,loJl-m{i)J-lo-l-m{-i) (5.8) 

I 

in order for Eqs. (|3.13|) and ( |5.71 ) to be equivalent. The preceding equation was verified 
numerically and the exact value of the constant is given in ( |4.32| ) . 

It is not possible to find an analytical expression for the coefficients /S.{qx = loG,m). 
Thus, one has to solve Eq. ( |5.1| ) numerically and then calculate every physical quantity as 
a function of the coefficients A{qx = loG,m). 

As in Sec. we construct the order parameter by taking a linear combination of 
the solutions of the hnearized gap equation. A{qx = (/q — 2p)G,m) has the form of a 
strip extended in the x direction and localized in the perpendicular direction around z = 
{—lo/2+p)c. Thus, the triangular vortex lattice with periodicity = Nc can be constructed 
as in Sec. 0, which leads to 

A(x, m) = A ^ 7pA^_^^e^('°-^^)^^ , (5.9) 
p ^ 
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where = A(g^ = lQG,m) is the normahzed solution of the hnearized gap equation (p. 11): 
J2m l^mP = 1- The amphtude A is chosen real. The coefficients 7p were introduced in 
Sec. The order parameter ( |5.9| ) can also be obtained from Eq. (41) of Ref.El by the 
appropriate gauge transformation and using A^ = A^^'"*^. Contrary to the QLA, the order 
parameter ( p.9[ ) correctly describes the entire phase diagram, from the weak field regime 
[ujc <^ T) where the superconducting state is an Abrikosov vortex lattice up to the very 
high field regime {uc ^ tz) where the superconducting state is a Josephson vortex lattice. 
The above form of the order parameter, obtained by taking a linear combination of the 
solutions of the linearized gap equation, is expected to be a very good approximation as 
long as the region of interest in the phase diagram is close to the transition line. 



B. Ginzburg-Landau expansion 

In order to derive the thermodynamics close to the transition line, we have to perform a 
GL expansion of the free energy. The situation is more complicated than in the QLA, because 
the supercurrents do not vanish and should therefore be taken into account. The Gor'kov 
equations have to be solved in presence of a non uniform flux density B(r) = H + b(r) where 
b(r) = V X a(r) is obtained from 47rj(r) = V x V x a(r). j(r) is the current distribution in 
the ordered phase. As noted by Rasolt and Tesanovicjl a consistent calculation of the free 
energy to order A^ requires the calculation of the Green's functions of the normal phase to 
ffist order in a(r). Such a procedure, although possible in principle, greatly complicates the 
calculation and will not be followed here. Fortunately, the situation gets simpler in the phase 
N = 2. In this phase, the supercurrents are of order t^/o;^ (Ref.0) and their contribution to 
the free energy is of order t\/uj'^. Thus, to ffist order in tl/ul, the screening of the external 
magnetic field can be ignored in the calculation of the free energy. This approximation is 
not justified in the phases iV > 4. We hope that it will give reliable results, at least not too 
far from the reentrant phase. Indeed, it will turn out that this approximation is sufficient to 
understand the main features of the superconducting phases. Nevertheless, we will see that a 
correct description of the phases iV > 4 would require a proper treatment of the screening of 
the field by the supercurrents. We could try to improve the crude approximation B = H hy 
assuming that the flux density B is uniform but not necessary equal to the external field H. 
This would however not lead to any significant improvements of the results. In particular, 
in the middle of each phase where the averaged flux density vanishes, this approximation 
becomes equivalent to completely neglecting the screening of the field. 

As in Sec. we exploit the translational symmetry of the order parameter using the 
magnetic Bloch states (Eq. ( |4.8| )). If we ignore the screening of the external field, the 
GL expansion can be written in the form (^4.20|) with 



A- 



S 





A.(r) 




A 



T ^ |A-(1,2)/Ap 



n T 1 Ag(l,2)Ag(l,3)*Ag(4,3)Ag(4,2)- 

where A|(r) = — Aj(r) is the variational order parameter (|5.9|) . Again, we use the notation 
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i = (qj, li). In (|5.10|) , the sums over are restricted to the MBZ defined in ([4.9| , |4.10| ). The 
pairing amphtudes A"(i,j) are defined by ( [4.13|) . Using (|5.9|) , we obtain 



J iV 

h—p^—m 



[at) . 



(5.11) 



As in the QLA, the total Bloch momentum of two paired states is equal to /qG. But the 
pairing is not diagonal any more in the branch index I, i.e. (f)°i is not coupled only to 
'PfoG-ci-io-i- Using the property (p integer) 



A^(g„ h + pN, h - pN) = e^^P^^'^^ A^(g„ h) 



(5.12) 



we obtain 



N 



a 



/3 



N 



A^(g., 



h,l2) 



A 



x((/o + h + l2)G, n) - x((/o + h + l2)G, T) 



where 



1 1 

E -r^A^(g.,/i,/2)A^(g.,/i,/3)^ 

h = lhMM 



T 



(5.13) 



-^4(^1: h) — — ^ 
ij 



1 



X ■ 



1 

iuj — avkx — al4^Uc){—iuj — avk^ + a(/o + l2)^c) 



(5.14) 



Unless otherwise specified, the sum over the integers li runs from —00 to 00. In ( |5.14| ), the 
sum over fc^. runs from —00 to 00. xi^xiT) is the pair susceptibility in zero magnetic field 
at the temperature T. Writing a = a'(T — Tc) with a' = da/dT\T^, we obtain, using ( p.5| ) 
and \l/(2;) ~ \n{z) - l/2z for \z\ > 1: 



a 



N{0) c Ag(g„/,-/-/o) 



We also have 



1=1 



N 



(5.15) 



/5-/3Bcsf EE 



1=1 



A-(g„Z,-/-/o) 



(5.16) 



The corrections to the above expressions of a' and /? are of the order of T^/(^c (^P 
logarithmic corrections for /3). In the quantum regime Uc ^ T, it is therefore sufficient to 
take into account only the primary gaps in order to calculate the coefficients a' and /3. The 
secondary gaps, which open at nuOc/2 {n 7^ 0) above and below the Fermi level (see Sec. 
|V (J| ) , do not contribute to a' and /3 to leading order in Tc/uc- However, they have to be 
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taken into account in the calculation of the critical temperature Tc- Comparing the primary 
gaps in Figs. [T3| - |T6t obtained numerically from ( ^.11 ), with those obtained analytically in 



the QLA (Eq. ( ^. 16| ) and Fig. |TT]), it turns out that the result obtained in the QLA is very 



close to the exact result. This means that the expression of the coefficients a' and (3 given by 
( [4.23| ) is very accurate. Thus, if we calculate the free energy and the specific heat jump at 



the transition, we will recover the results of Sec. the only difference being that the value 
of Tc is now exact. Again, we find that the minimum of the free energy in the reentrant 
phase is minimized for N = 2. When the field decreases from its value in the reentrant 
phase, the system undergoes a first order transition and the minimum of the free energy is 
then obtained for = 4. As in the QLA, the best value of N switches to 6 before reaching 
the next first order phase transition. This result, which seems to contradict the assumption 
that the first order phase transitions are due to commensurability effects between the period 
of the order parameter and the crystalline lattice spacing, indicates that the contribution of 
the supercurrents to the free energy has to be taken into account in the phases > 4. 

C. Excitation spectrum 

The excitation spectrum can be obtained from the BdG equations.0 Performing the BdG 
transformation 

C(r) = E(<(r)7:,. - cx<(r)*7:y , (5.17) 

n 

and expanding the coefficients u and v in the basis of the magnetic Bloch states: 

<(r) = E<.C(r), 

i 

<(r) = E<.C(r)% (5.18) 

i 

we obtain the following BdG equations: 

I' 

iK,a + e^G-q,/,^)<.oG-q,. - E <^,l'^Ulz. l'^ 0* = . (5.19) 

I' 

In the quantum regime, uJc ^ Tc so that ^ A even at T = 0. The latter inequality allows 
us to treat the pairing amplitudes A"(g2, 1, 1') perturbatively. At low temperature, the order 
parameter ( |5.9| ) is modified by contributions of superconducting condensation channels with 
critical temperature < Tc. However, the main contribution to the order parameter still 
comes from the channel with the highest critical temperature so that ( |5.9| ) should remain 
a good approximation. To leading order in A'^/Uc, the effect of A"(g^,/,/') is to lift the 
degeneracy between and e^Q_q ;/ g:. Thus, up to corrections of the order of A^/tu^, a 
gap 2|A"(g2, /, /')| opens in the spectrum at qx = akp — {I — Iq — l')G/2 (with the restriction 
that qx has to be in the MBZ). This gap is located at a(/o + / + l')u}cl2 away from the Fermi 
level. The primary gaps A"(g2, /, —Iq — I) (the only ones which were considered in the QLA) 
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open at the Fermi level, while the secondary gaps open above and below the Fermi level. 
The resulting spectrum, shown schematically in Fig. 17, is very reminiscent of the one of the 
FISDW phases.BaEJ Since the one-particle states are localized in the z direction on a length 
of the order of ct, the pairing amplitudes A"(g2, /, /') are important only if |/ — /'| < t. This 
means that there are ~ ~ t secondary gaps with a significant value opening above and 
below the Fermi level and extending on an energy width of the order of tz- In the reentrant 
phase, the secondary gaps are of the order of with respect to the primary gaps. 

When the field decreases (at T = 0) below the semiclassical critical field ifc2(0), the 
amplitude of the order parameter will grow so that A will become larger than 10^- The 
coherence length f/A then becomes smaller than the (longitudinal) magnetic length 27r/G 
and the quantum effect of the field (i.e., the bending of the semiclassical orbits by the 
field) can be ignored. Thus, the condition A ~ (at T = 0) signals the crossover to the 
(anisotropic) Abrikosov vortex lattice state. 

As noted above, the QLA very accurately describes the primary gaps and therefore 
provides a very good approximation of the minimum excitation energy. We also note that 
the zeros of the primary gaps which were obtained in the QLA are not destroyed by the off- 
diagonal pairings (/, /') (/' 7^ — /q — /). A similar result was obtained in the case of isotropic 
superconductors where it has been shown that the gapless behavior obtained in the QLA is 
not destroyed by a weak off-diagonal Landau level pairing.1'0 

As in the QLA, we obtain a non monotonous behavior of the minimum excitation energy: 
it decreases with the field and vanishes in the phase = 8 at = L7 T, but becomes 
finite again in the phase = 10. This gapless behavior turns out to strongly depend on the 
dispersion law of the non-interacting system and is therefore accidental. For example, if we 
add a second neighbor hopping term i!^ co^(2kzC) to the dispersion law (|2.1|) , the zeros in the 
excitation spectrum appear for different values of the magnetic field. A proper treatment 
of the screening of the magnetic field would modify the BdG equations and is expected 
to suppress this accidental gapless behavior. Nevertheless, when the field decreases, the 
minimum excitation energy decreases and becomes very small in the large A^ phases as can 
be seen in Fig. 18. This figure shows the gap which opens at the Fermi level in the phases 
A^ = 26 and A^ = 28 (_the order parameter A^q ^^(1") corresponding to the phase A^ = 26 is 
shown in Fig. 7 of Ref.Ej). Even if we believe that the gapless behavior of the phase A^ = 26 
is accidental. Fig. |1^ clearly shows that the large A^ superconducting phases effectively 
become gapless (note that one can distinguish 21 different phases (i.e. N = 2,- ■ -,42) in 
the quantum regime: see Fig. lb of Ref.lli). It is however difficult to conclude if the system 
evolves towards a real gapless behavior. Let us finally mention that the wrong estimation 
of the minimum excitation energy in the phases A^ > 4 is related to the poor estimation of 
the free energy in these phases. 



D. Current distribution 

In this section, we give the expression of the current distribution in the superconducting 
phases. We will show that the primary gaps do not contribute to the super currents. Since 
the calculation is very fastidious, we only give the final expressions. The current distribution 
is obtained from 



28 



j^(x,m,m + 1) 



ect 



2i 



(5.20) 



Here, jz{x, m, m + 1) is the current at point x between the chains m and m + 1. is the 
correction to the Green's function which results from a non zero order parameter. To lowest 
order, 



SG^ir, r', iu) = -Y: 0?(r)0^(r')*A^(l, 3) A^(3, 2)*G^(1, u;)G§(3, -a;)G°(2, u) , (5.21) 

1,2,3 

where 0f is the magnetic Bloch state defined by ([4.8|) and we use the notation i = (qi,/j). 
Note that to lowest order (A^), the screening of the external field can be ignored and the 
correction can be calculated with the Green's functions G"(i,co') of the normal phase 
in the presence of a uniform magnetic field H. From ( |5.20| , |5?21| ), we obtain the following 
expressions for the Fourier transform of j(r) at wave vector k = {piNG,P22tt/Nc): 



N 



jx(,Pl,P2) 
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+ (g„/i,/3)A+(g„/2,/3)*, (5.22) 



where the function F/v m is defined by 
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if AT = M , 
if AT ^ M . 



(5.23) 



The function F^^m has been introduced in Ref.El where its expression is given by Eq. (50) 
for ujc ^ T. If we retain only the primary gaps A"(g2,Z, — / — Iq), then F = Fq^ = so 
that the supercurrents vanish. Thus, in the QLA where the secondary gaps are ignored, 
the supercurrents vanish. The current distribution can also be obtained as a function of the 
coefficients A^ using the relation 



1 



N 



P}J{lzC-p^){h-l2) 



xe-*5('i+'2-2Hj^^^,^_2^|^2tsin 
Again, we only give the final expressions: 



(5.24) 
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^/ g-i«(pi7V-/i-2mi+Z2+2m2)j / _2tsin 



U+P2 



N 



J/2 (-2t sin m) )r7pj,p2 , (5.25) 
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N 



\Ji^{-2ismu))r]p^^p^, (5.26) 



where we use the notation (■•■) = /o'' ■ ■ -f^- The function ?7pi,p2 is defined by 

{1 if pi and p2 are even, 
i if pi and p2 are odd, 
otherwise. 



(5.27) 



Since Aj^ = A^^'"*^, Eqs. ( ^I25| , ^I26D are analogous to Eqs. (48,49) of Ref. Ill where the current 
distribution has been calculated in the gauge A'(0, 0, —Hx). The current distribution was 
calculated numerically in Refill using ( |5.25| , O^ ). In the reentrant phase, ( |5.25| , p^2^ ) can be 
simplified by retaining only the terms of order t^: the current distribution is characteristic of 
a triangular Josephson vortex lattice. In the phases iV > 4, the current distribution shows 
a symmetry of a laminar type and is different from what is obtained in the Abrikosov or 
Josephson vortex lattice.El 



VI. CONCLUSION 

We have presented a systematic study of the phase diagram of a quasi-lD supercon- 
ductor in a high magnetic field. We have obtained the thermodynamics quantities and the 
quasi-particle excitation spectrum of the quantum superconducting phases which are sta- 
bilized at high magnetic field as a result of the magnetic-field-induced confinement of the 
electrons. The reentrant phase (very high field limit) is the natural limit for the study of 
superconductivity, since the orbital frustration of the order parameter is suppressed and the 
screening of the external magnetic field can be ignored, which considerably simplifies the 
analysis. Although we have not taken into account the screening of the field in the other 
phases, our analysis clearly shows how the properties of the system evolves when the field 
is decreased from its value in the reentrant phase. 

The main results can be summarized as follows, i) In the reentrant phase {uj^ ^ tz), the 
superconducting state is a Josephson vortex lattice. The behavior of the system is very close 
to the zero-field BCS situation, up to corrections of the order of tl/ojl. For instance, the 
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supercurrents arc of the order of t\luj\ and the specific heat jump at the superconducting 
transition ^C/Cn is very close to the zero-field BCS value (AC/CAr)Bcs- The phase is 
paramagnetic due to the positive slope of transition temperature {dTc/dH > 0). Gaps 
open at the Fermi level on the whole MBZ. The dispersion of these gaps is of the order of 
tl/cu'^ and the minimum excitation energy is finite. Gaps also open at ncUc/'2 {n integer) 
away from the Fermi level but are of the order of tl/ojl with respect to the gaps opening at 
the Fermi level, ii) In the other phases {N > 4), the behavior of the system deviates more 
substantially from the zero-field BCS situation. The phase is first paramagnetic and then 
diamagnetic, which is a consequence of the sign change of the slope of the transition line 
Tc. The specific heat jump at the transition becomes smaller than the zero-field BCS value 
{AC/C]^)bcs- There are ~ iV ~ t gaps with a significant value opening below and above 
the Fermi level. The dispersion of the gaps which open at the Fermi level increases. The 
minimum excitation energy decreases and the quasi-particle spectrum becomes gapless for 
large N. 

Our result have been obtained within a simple model where the electrons interact through 
an effective local attractive interaction. Some of these results would be modified in a more 
complicated (or realistic) model. For example, in the case of a rf-wave (with respect to the 
X and y axis) superconductor, the excitation spectrum would be gapless for any value of the 
field. The inclusion in our analysis of the PPB effect, which yields to the formation of a 
LOFF state, would also result in a gapless behavior. However, the general structure of the 
phase diagram does not depend on the details of the microscopic mechanism which is at the 
origin of the attractive electron-electron interaction: the existence of superconductivity at 
high magnetic field in a quasi- ID superconductor originates in the magnetic-field-induced 
dimensional crossover and does not rely on a particular model of superconducting pairing. 
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FIGURES 



FIG. 1. Diagrammatic representation of the ladder approximation for the two-particle vertex 
function r°"'(ri, r2; r'^^, rg). The zigzag Hne denotes the attractive electron-electron interaction, a 
and a' refer to the sheet of the Fermi surface. 

FIG. 2. Solid lines: critical temperature vs magnetic field for Iq = and Iq = 1 in the QLA. 
Dashed lines: exact mean-field critical temperature (see Fig. 1 of Ref. 19). Tco = 1.5 K and tz = 20 
K. 

FIG. 3. a) Self-energy correction in the Born approximation due to impurity scattering, b) 
Vertex correction due to impurity scattering for the pair propagator. The dashed lines with a cross 
denote impurity scattering. 

FIG. 4. Amplitude and phase of the order parameter A/Q^^(r) (Eq. (|4.2|)) obtained in the 
QLA in the phase Iq = 0, N = 2 {H = 5.8 T). In order to compare with the exact order parameter 
obtained in Ref. 19, we have made a gauge transformation from A(i/z,0,0) to A'(0, 0, —Hx). 

FIG. 5. As in Fig. |, but for the phase Iq = I, N = A {H = AT). 

FIG. 6. As in Fig. g but for the phase Iq = {)^N = Q{H = 2.4 T). 

FIG. 7. As in Fig. § but for the phase Iq = I^N = ^{H = 1.7 T). 

FIG. 8. As in Fig. |, but for the phase Iq = {), N = l{) [H = 1.3 T). 

FIG. 9. Excitation spectrum ; ^ = / of the normal phase in the representation of the 
magnetic Bloch states i;^^;. qx is restricted to the MBZ 
N = 2. 

FIG. 10. Ratio r = (AC7/C^)/(AC/C7v)bcs vs magnetic field. 

FIG. 11. Pairing amplitudes |A°(g2,/)| in the QLA in the phases AT = 2, iV = 4, iV = 6, 
N = 8 and N = 10. The units are chosen so that max |A°(g2, ^)| = 1. 

FIG. 12. Quasi-particle excitation spectrum E^^^ in the QLA in the reentrant phase Iq = 0, 
N = 2 for the branches / = (a) and I = 1 (b). q is restricted to the MBZ. The gap opens at the 
Fermi level. 
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FIG. 13. Pairing amplitudes | A°(g2, 1, L — I — in the exact mean-field analysis in the phase 
N = 2 {H = 5.1 T). A gap 2| A"(g2, l,L — I — opens in the branch I at aLoJc/2 away from the 
Fermi level. The units are chosen so that max | A"(g2, l,L — I — Iq)\ = 1. 



FIG. 14. As in Fig. y, but for the phase iV = 4 (i? = 2.8 T). 



FIG. 15. As in Fig. 13, but for the phase N = Q {H = 2T). 



FIG. 16. As in Fig. 13, but for the phase AT = 8 (i? = 1.5 T). 



FIG. 17. Schematic representation of the quasi-particle excitation spectrum for the right sheet 
of the Fermi surface, qx is restricted to the MBZ. Gaps open at nuJc/2 (n integer) away from the 
Fermi level. The dispersion with respect to qz is not shown. 



FIG. 18. Pairing amplitudes (primary gaps) |A°(g2, /, — / — /o)| in tlie exact mean-field analysis 
in the phases N = 2Q{H = 0.48 T) and N = 28 {H = 0.45 T). 
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